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ABSTRACT 


It has become evident in seismology and seismic 
prospecting that to obtain more consistent results with 
actual field data a more sophisticated mathematical model 
than that of an isotropic homogeneous earth must be 
considered. To this end asymptotic ray theory applied 
to anisotropic inhomogeneous media, has been presented 
by Cerveny (1), with the special case of transverse 
isotropy dealt with in the second of his quoted papers (2). 
With the object in mind of using these results to produce 
synthetic seismograms for transversely isotropic plane 
layered media it becomes necessary to compute first, 
reflection and transmission coefficients of the coupled 
quasi-compressional (P) and quasi-shear (SV) waves at an 
interface of two such media. These aearetcientss along 
with the free interface case will be presented here. 

As the requirements for continuity of displace- 
ments and stresses will be developed in terms of asympto- 
tic ray theory, a brief treatment of it (after Hron and 
Kanasewich (6)) will be presented first. Elements of the 
theory of Cerveny's papers (1),(2) will be quoted, in so 
far as it is necessary to obtain expressions for the 


normal velocities of the wavefronts. 
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CHAPTER 1 
ASYMPTOTIC RAY THEORY FOR AN ANISOTROPIC 


INHOMOGENEOUS MEDIA 


In general the equations of motion in an elastic 


media, in Cartesian coordinates are (neglecting body 


forces) 
0; 
iar) pu, = a, 


where Aes uot is the internal stress force, op is density 
and u, is acceleration. 


The stress tensor O55 is related to the elastic modulus 


tensor Chie (which is in general coordinate dependent) 
as 
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where 


satisfy the conditions of 


The elastic parameters Ci akg 


symmetry 
ee ak ee ig 
reducing the number of independent elements of the elas- 


cic modulus tensor from glato 21. Supstacucing tl. 2) 
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into (1.1) and using the symmetry properties of the 
Ci sk yields 
(1.3) | ae bee 
5 OX; IK Le x 
It will be assumed that the solution of (1.3) 


near the wavefront can be written generally as 
(1,4) u, (x,/t) = W, (x,t) s(t-1(x,)) +... Wipe sda 28 


where s describes the source as a function of time and 
may be discontinuous at t = 0. W, is assumed analytic 
in the neighbourhood of the wavefront. The equation 


describing the position of the wavefront at time t is 
CoS ert) Th LZ 


tT being the phase function which depends upon the ray 
path from the source to the observer at M(x). 

Let At = t- T be the time a short distance from 
the wavefront and expand W; in a Taylor series about 
t = t+ At. This Taylor series expansion is justified 
only if there are no fluctuations of great magnitude in 
the elastic parameters of the medium in a distance of 
one wavelength.” As° it» will be shown later that the 
normal wale chley to the wavefront is dependent on the 
elastic parameters, the expansion can be assumed valid 
if the gradient of the normal velocity is slowly 


varying. (1.4) becomes 
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{ } is totally spacially dependent while [ ] is time 
frequency and spacially dependent. 


A monochromatic sinusoidal solution of (1.3) is 


Cryo) u; (x rt) = Be eeees? Mee 


Assuming a; can be expanded in a power series in iw 


yields 
fore) ia 
Bene ) n 
n=0 | Cru) 
and (1.6) becomes 
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A particular case for harmonic waves is when s in 


equation (1.5) is 


n iw(t-T) 
ate s(t-1t) = = arihg are) 
i (iw) 
Equation (1.6) is in a useful form since the at are 


independent of frequency with the frequency dependent 
position of both the source function and wave modifica- 
tion due to the medium contained in the term Ss, (t-t(x,)). 
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Equation (1.8) can be generalized to an impulsive source 
function by first obtaining the Fourier transform of the 


desired impulse response s(t) 
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i.e. S(w) = | s(tye tt ae . 
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After summing (1.8) ovér all frequencies, the 


displacement at any location is given by 


u, (x,rt) = AE ) ite | = 
n=0 if (1W 
O 


~ 


as S(W) is required to be negligible for w < oe (see 
high frequency approximation when the low frequency part 
of the source spectrum is assumed to be small compared 
with the rest of the spectrum). It should be noted that 

ds_ (&) 

n 
Nt ieee Sy-1 (5) : 
Subetleucing. (1.6)) into, Cl. 3) cyields upon mani- 


pulation of the summations, along with the additional 


assumption that ee =a,7t) = OF t= ha, 3, that 
ie Goh ee (n) 
(1.9) ett) ta, 255K 2PiPeAK } 
(n=) 
oA 
k dy fone) (n-1) 
+ tas 5K Pa ox, 5 Ox, (Pas 5p pP pry 
ae 
Sire (pa jk} = 0 
p Ox. ijk OX p 


where Chats Sign?” and Pp; = dT/ox, the slowness vector. 
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The p,'s at a given point and time are related 
to the normal velocity of wavefront at a given point 


and time by the relation 


N, being the directional cosines of the wavefront normal 
and. Vv the nermal velocity. “For (1.9) to be satisfied, 
the square bracketed term must be equal to zero for all 
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By definition the tensor a. has the symmetry 


LIK 
property 


Panes | KR a" 


Therefore Tak is symmetric and as a result the eigen- 


values of Dak are all real. Since the strain energy is 


positive definite, Tat is positive definite and thus 


the eigenvalues of Tsk 


From this it can be seen that 


are all positive. 


fel oil 2a) det(lay - G6 08 


sb 
yields 3 real and positive eigenvalues Gis Gor G3. One 
eigenvalue corresponds to a quaSi-compressional (P) 

wave and the other two correspond to two independent 
quasi~shear waves, SV and SH. 

The only non-trivial solution of a 8 53) Ago = 0 
occurs when one of the eigenvalues Gor sie 1, Zyo- (On Dt 
is equal to unity. Neglecting degenerate cases, there 
are three possibilities, corresponding to the propaga- 
tion of one of the three different types of wavefronts. 
It is from these eigenvalues using the relation pai 
Ni/Vor s = 1,2,3, that the normal velocities to the 


wavefronts are found. 
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CHAPTER 2 


THE SPECIAL CASE OF THE TRANSVERSELY ISOTROPIC MEDIUM 


Because of the symmetries of the a. a more 


Ee 
convenient notation arrangement may be adopted. Asa 


result 255K Amn according to the following scheme: 


The elastic parameters of a medium can then be tabulated 
in a 6x6 symmetric matrix with 21 independent entries. 
However, any calculations with all of these 21 parameters 
becomes unwieldy so some simplifying assumptions must be 
made. Cerveny (2) suggests a model of an isotropic 
inhomogeneous medium which can be described by the 9 


independent parameters. 


eee) | en ee J : 
al Gye mth) bl ° 4 
Pig mie 300 eos ‘ 
0 0 @ Maye 0 0 
0 0 0 Nese yo 1 0 
0 0 0 0 6 PaRS 6 


é {33 SMS Aas ee bo el. i at oe 


wb stuont, ah 


t? hy je 
sii @ ete a edt oe ens 20 
“ 4 oan 
n hoe 
* ee 4) Ye 
; 4 j : 
; f 
" PS ie \ 
a crue & 
heteiodes od nena nav jap chomaagy ache & 
. aid : 5 det JmoeLrag< hw ris rs atie ti 


rT _ Tun, 


iweb 7 
i ; { 7 Aw 
ot yen Brio he 3. OR i 2 ‘atin ie A rahe Be 
{ an estichatimaie as. “9s a Teint ag oe eh _ 
St SRA Ome ee Ao. Labames ataeirs a3 (S) siesta 
y ae ya bedi ieee 4a); ris oi pane deeb 


A special case of the model, transverse isotropy, 


will be considered here, for which A =A 


lt Seep eee 55° 


Ay 3 = Aj 3 Ajo = Ali 2A66- Inherent to this special 

case are the additional assumptions: 

{a) The elastic parameters depend on only two 
coordinates which specify epicentral distance 
(x), sand, depth .@je 

(b) The slowness vector has only two components, Po 
chosen to be identically zero. 

Under these assumptions, the source lies in the plane 

(x,Z) as does the ray under consideration for all time 

aS well as the corresponding normal to the wavefront. 
For transverse isotropy the matrix elements Psy 


are given by 
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The solution of (1.12) yields the three roots 
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G, = 1 corresponds to a quasi-compressional (P) wave- 
front propagating in the medium, while G. = 1 corresponds 
to a quasi-shear (SV) wavefront and G, = 1 corresponds 


to a quasi-shear (SH) wavefront. 


The normal velocities Vo and Voy of the P and SV 


type wavefronts can be obtained from (2.2) with Go = l, 
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for the quasi-compressional case and 


(2.4) Sin 9 ce COS. 
SV SV 


for the quasi-~shear case, 9 being the angle between the 
normal to the particular wavefront at some time and the 


Z axis. The normal velocities are given by 
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the eigenvalue problem are 
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CHAPTER 3 


REFLECTION AND TRANSMISSION COEFFICIENTS AT AN INTERFACE 


OF TWO TRANSVERSELY ISOTROPIC MEDIA 


3.1 Ineroduction 


In this section a general set of linear equations 
using asymptotic ray theory to calculate reflection and 
transmission coefficients in the case of a curved inter- 
face, between two transversely isotropic inhomogeneous 
media in welded contact will be developed. 

Let ) be an invedeacs Ge first order, and let 
be analytic. Denote the point at which an arbitrary 
wavefront strikes the interface ) by 0. The plane of 
incidence will be defined as the plane defined by the 
tangent to the wavefront normal and by the normal to the 
interface at the point 0 (see Fig. Al). It is assumed 
that the incident wavefront is propagating in the upper 


medium (medium 1). 


Some Notation: 
The displacement vectors arriving at or leaving 


a boundary are specified by 
> > 
uj (x,y,2,t) = a haay Guy aye, (tar) 


where the index v specifies the phase and medium 


Lz 


poe 


sno heups Resad fh. ie tee: kSReueDp 8 

hang nabisigifitios 

~ sapees beyaw> 6 
Riloomeponortni | 


boyolavah. ect fdiw eet benitew ih ‘ 

< Bal baw .2Shxo Jer ae ro Parise list a4, ad J tea 
visit teas m auitw Js der dqc; = 4 Er * ggon80 
to onsiq aa? »?! Ye % portage ti “iid ohio 


eile ud eink? ob Sela alts! ag peal tap sel Me 


7Sqqu oct mi coivepsgeig ei 1 ohcialecalaaal: 


OaivEod to ds prise avoseey an 


incident plane 


normal to 
interface 


tangent to 
wavefront normal 
at O 


incident 
wavefront 
normal 


interface 2% 


Fig. Al Geometry of incidence. 
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: : : gh () 60) korg > 
=0 incident wave (either P or SV) A =a ors 
re) no no 
=l reflected P wave xn) _ p 
at nl 
(3.1) « v=2 transmitted P wave Ain) Ie 
=3 reflected SV wave x(n) 2 p 
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In general, for a transversely isotropic medium, 
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Pate ot /Ox, Po = 0, P= 91/dz must be continuous at the 
interface. From the definition of Py the following 


condition results 


the v being defined as in (3.1) and the ce are the acute 
angles between the wavefront normal and the Zz axis. 
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approach zero as the wavefront moves away from the 


interface. Thus the minus sign is chosen. 


3.2 Incident Quasi-compressional P Wavefront 


As shown in Fig. A2, the solid lines indicate nor- 
mals to the P-type wavefronts, while the dashed lines de- 
fine SV type wavefront normals. This notation will be 
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used throughout the thesis. The eigenvectors ty are 
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Fig. A2 Orientation of wavefront normals and 
eigenvectors for an incident quasi-compressional 
P or quasi-shear SV wavefront. 
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The conditions which must hold at the interface 
between two layers in welded contact and the continuity 
of the x and z components of displacement, and the con- 
tinuity of normal and shear stresses. Introducing 
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3.3 Incident Quasi-shear SV Wavefront 
The geometry of incidence is shown in Fig. A2 
and the four continuity equations obtained in an analo- 


gous manner to the preceding are: 
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where Vo denotes the normal velocity of the incident SV 
wavefront. 

The eqwectons -C3.i/)r7 "GS. 6)", (3..9).).(3..L0) and 
C3 elle); (3.12), (3.13), (3.14) are two recursive sets of 
four linear equations in four unknowns. To know the so- 
lution of these two sets for a given n, the solutions for 
all m<n must be known. 

If the geometrical parameters of the layers are 
much greater. than the wavelength, and if the observer is 
not very close to the source of the wave, taking only the 
first term in the asymptotic series (n=0) can be justi- 
fied to give a reasonable approximation. This so called 
zero order or plane wave approximation will be considered. 

As a further simplification, it will be assumed 
that the axes of anisotropy in both media are perpendicu- 
lar to the interface. Thus, for an incident P type wave- 


Eront 0. F @ and for an incident SV type wavefront 


ik y 
Po 83. 


23 


we YS 

SS taal 
A ; , 
' vii } ’ 
78 u 

aa ‘n 


hi a afe Saaed Fi 
: 


U2 anehioad ont to vilop Tew Losin 


bas , (er. BY’ she .€ vel is wal Sea i 

Nick ater’ midiexvonx wit ove (KEPT (i ae £) Ags BY) (2t8), aan 

| Oe erid woud igre’ ery oe sae nt Lshwas iad tl aneriss | a 

zor -pnoiau¥oe i a hipints B ‘tot eta ows sands to r 

a od Jaume in ie, 

O28 ieee a itd to otktuadnd ihietbied one ar, | 7. os 

ek sevredtic ond ti Sms sty ae leven ad? cnt asteexp baie > 
oct yino patdet , avaw, orld’ to eax@od oe ‘oa agate’ grav gon 

“itept ‘ou fan cee charlie ‘alsomeiyps ert nk ore jer, a 

bellso o@ eit .aadvsmd sq | sidanowess & avap a gis 8 

Bessbieno eb LE bw soRpenaomas avew reese to sebze 

irae ie ah tk ee. WRN, ur 


Let sin 6, = X. All other pertinent angles will be para- 


meterized in this variable via Snell's Law. 
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Using these above formulae, the following simpli- 
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After implementing the substitutions in the zero 
order equations, the two following sets of four linear 
equations in four unknowns result. The unknowns are 
the amplitudes of the reflected and refracted disturb- 
ances, as the amplitude of the incident disturbance is 


assumed known. 
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Incident P type wavefront 
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Selvine’(3.15) and (3.16) using Cramerts 
method with the incident displacement set to unity the 
rerlection and transmission coerricient Ay tabulated 
below are obtained. (According to our convention, 


uw = 1 corresponds to an incident P wavefront and 


U 5 corresponds’ to an incident’ SV *wavetront., while 
the second subscript v typifies the resultant ref- 


lected or transmitted wavefront.) 
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3.4 Free Interface 

The free interface case has a wavefront (either 
quasi-shear SV or quasi-compressional P) propagating in 
a transversely isotropic medium and impinging on the 
interface of that medium and a vacuum. Reflected shear 
SV and compressional waves only result. 

The boundary requirements are the continuity of 
shear and normal stresses along with the previously dis- 


cussed conditions on the first derivatives of ee (see 


FigeiA3)% 
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The requirement for the continuity of shear stress 

(a) + (a5) + (i,) = 0 yields the following two 
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Fig. A3 Wavefront normals and eigenvectors for 
the free interface case. 
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Equations (3.18) , (s.19)) and (3.20);7; (3.21) are two 
sets of two recursive linear equations in two unknowns 
whose solutions give the required reflection coefficients. 

In the zero order approximation (again assuming the 
axis of anisotropy is perpendicular to the interface) these 


systems become (3.22) and (3.23) after the following sim- 


plifying substitutions. 
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Incident SV wavefront 
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efficients of reflection. 
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3.5 Surface Conversion Coefficients 

Apreceiver situated yom the earun'’s surface 1.¢€. at 
the interface of a transversely isotropic medium and a 
vacuum, aecords nov only the=disturbance* caused by “an 
incident wavefront at that point, but also the disturbances 


resulting from the two reflected wavefronts. 
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bete te Vdenote \the-=inverface ; iG ’thevlocavion vor 
the receiver at the interfacé and a” (G) the displacement 
registered by the receiver. (v = 0 or 6 depending on 
wnevner. the incident ..waveitront*vs of 'the- -P sor, SV type.) 


aus 
aay Magis? > > 
(3.25) Us a Gy eee) a Ce) 


where the first term on the R.H.S. of (3.25) is the 
contribution from the incident wavefront and the second 
and third terms. are contributions from the reflected P 
and SV wavefronts respectively. 

In the ‘zero order approximation in asymptotic 
ray theory an arbitrary displacement vector can be defined 


by 

wrecG) = A (G)exp iw(t-1 (G))n (G) 

=} 

n,G@) being a unit vector. 

As T UGy) = T,(G) = Tt (G) at the interface, for an incident 
P wavefront (3.25) becomes 


(3.26) a” (G) = Poo(G) exp iw(t-1,(G)) {m5 (G) 
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and similarly since Ta (Ge = T 5 (G) = T,(G), the resultant 


expression for an incident SV wavefront is 


(3.27) wu” (4) = pee, exp iw (t-tx(G)) n= (G) 
" Re ot, (G) + Ren), (G)} 


The R are given by (3.24). 
> ee POT > 
Let g,(G) ee) + Ron, (G) + Ryn, (G). 
g(a) is called the surface conversion vector and the x 
> 
and 2z% components of B(s), Bey 
surface conversion coefficients. From (3.24) it follows 
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CHAPTER 4 


ELLIPSOIDAL VELOCITY ANISOTROPY 


4.1 Introduction 


Analytic expressions for reflection and trans- 
mission coefficients at the interface of two traversely 
isotropic media were obtained in Chapter 3 using the 
zero-order approximation of asymptotic ray theory, which 
amounts to using a plane wave approximation. Since the 
expressions involving the wavefront normal velocities 
are still quite complicated, and as these velocities 
are required in calculations involving Snell's Law, a 
Simpler model is sought to analyze here. The model which 
will be considered is one consisting of homogeneous plane 
layers which display ellipsoidal velocity anisotropy (all 
wavefronts are ellipsoids of revolution). The major axes 
of these ellipsoids will be assumed to be either parallel 
or perpendicular to the plane interfaces. 

As was shown in Chapter 3 a transversely isotropic 
medium can be specified by the five modified elastic 
parameters Ay. A, 3: A337 Acer and Ace where Ai; ad Cis /p" 
0 being the density of the medium, which is constant if 
the medium is homogeneous and C3 are the elastic para- 
meters of the medium. Since Ace is inherent to SH type 


wavefronts, which are not considered, only 4 modified 
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elastic parameters need be known for the purposes of 
this discussion to describe the medium, The further 
constraint that the wavefronts are all ellipsoids of 
revolution reduces this number to 3, since for ellip- 
soidal wavefronts the condition 


2 


55) 7 'Ay, 7 Ags) (A337 Ags) = 0 


(42) (A, 3+ A 


holds, and as a result A 3 can be expressed in terms of 


ei 


A and A (see Gassmann (5)). 


Lae 33! 55 
Upon esubstituticon of (4.1) anto (2.4) and (2.5), 
the equations for the normal velocities, the following 


expressions are obtained. 


2 


| ag n nae 
(4 42) V. ™ A. sin 6..+-A.,cos.6-= A,3+ (A,,7~ A33)sin fs) 


P 1 a3 


(4,3) V = VA 


4.2 Ray and Normal Velocities 


The angle and velocity that are most often con- 
sidered in seismological applications are the ones 
describing the ray and not those of the wavefront normal 
associated with the ray. “It is the latter of these two 
which has been considered until now, and except for 
spherical wavefronts, the two are in general not the 
same. As can be seen from equation (4.3) the normal 


velocity of the SV type wavefront is independent of any 
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angle, from which it follows that the SV type wavefront 
is forced to be spherical and hence the ray velocity 
and normal velocity are the same in magnitude and 
direction. 

Using Euler's theorem on homogeneous functions, 
the nonlinear partial differential equation G (Py als 
(m=SV or P) can be solved by means of characteristics 
(Courant and Hilbert (4)). The equations of the charac- 
teristics corresponding to the partial differential 
equation G (Py )=1 are given by 
(4.4) ioe 1 On Pil 

dt op; 
The quantities dx, /dt are the components of the ray 
velocity. 

With the simplification (4.1) these velocity 


components are 


Ax, gloat Ady Sin 6 

dt Kapa st V, 
(4559 

dx, Va A233 cos 6 

at 3 38 V, 


for P type rays, and 


J 


“ep breleaqomsedil 1 tne yd bev Ye neo (x 


eke atey? au Cae 


sted og? 2? nor 


(et add Fo) patria pea ety o16 nine sotsi2ae 


witsiinlts eset (1.b) noitgoihtignula exit oan 


4 anh He 


ae: r 


\ on ~) 
xn. fr ae 
7 < ¥ 


ne ‘AS 


i ae 


Se 


) 
° 


j > eae eae 


PS, @ ; 


d nage coe) 
Sr ne 


inh 


, . ; 
iy i 
n 4 Mi L i a 
i me. ps a 
utd Pe 


eu Btn 7 Ars sin 60 
dt So ie 
(4.6) 
ie oe aa, Ane Sin 6 
ae y 
SV 


for SV type rays. The angle 6 is again the angle the 
wavefront normal makes with the vertical (x5) axis. 

The ratio of the ray velocity components yields 
the tangent of the angle, $9, the ray makes with the 
vertical axis. 


A 


(4.7) tan ¢ = — tan 0 for P type rays, and 
33 
(4.8) tah ¢ = tan 06 for SV type rays. 


From (4.6) and (4.8) it is easily seen that the 
SV type ray velocity and angle are identical to the SV 


i pe normal velocity and angle. Using (4.5) and (4.7) 


ax ax 1 
tie P type ray velocity c = teen) + EEE aN can be 
at at 
determined. 
2 
a yy sin“ 4 cos*¢ a Vo 
poe 07 7 Sit ae aS sin O+A2 pees 
. LA 35 


Vp being the normal velocity of the P type wavefront. 
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The relation between wavefront normal angles and 
velocities and ray angles and velocities can be seen 
graphically in Fig. A4. 

It should be noted that at 6 = ¢ = 0 and 6 = $6 = 
1/2 the normal velocity and ray velocity are equal; i.e. 
co = We = YA33 = ag “anid c™/? |e Nay? = YAj1 = ax. It is 
assumed that these two velocities can be measured as 
can VAge = b, the SV ray and wavefront velocity, along 
with the density. With these parameters known, the 


medium can, for the puposes of this thesis, be consid- 


ered fully described. 


4.3 Snell's Law 


The simpler expressions obtained for normal velo- 


cities can be substituted into the statement of Snell's 


law 
sin 0 sin 6 
(4.10) pc tO we % 
v5 Vo 
v=0 - incident wavefront 


v=l - reflected P wavefront 


v=2 - transmitted P wavefront 


v=3 reflected SV wavefront 


v=4 - transmitted SV wavefront 


and if the angle of the incident wavefront is known, 


expressions can easily be obtained for the angle of the 
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reflected and transmitted wavefront normals. 

AS an example, consider a P type wavefront inci- 
dent from medium 1 on an interface. The angle of the 
transmitted P type wavefront normal in medium 2 is given 
by 


(4.11) sin 0, =az2 sin @./{azl* + (axl? - Pes eg ax2° + az2°) 
e 
x sin“6 ee 
fe) 
In a Similar manner all other wavefront normal angles 


can be expressed in terms of the sine of the angle of 


the incidant and/or reflected P wavefront normal. 


4.4 Numerical Results 


A program was written using Fortran IV which 
calculates reflection and transmission coefficients, 
including the free interface case and surface conver- 
sion coefficients, using the formulas developed earlier 
in the paper for varying angles of ray incidence. 

A Calcomp plotter was used to display the results 
which consist of two graphs for each coefficients; one 
of m@dulus of the coefficient vs. ray angle of incidence 
and the other with phase of the coefficient vs. ray 
angle of incidence, since in general the coefficients 
are complex. (The phase in the following plots may in 
some cases vary by 7, since in writing the program, the 


Z-axis was shown to point downwards.) 
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Since several situations are to be investi- 
gated space dictates that not all coefficients can 
be presented. The four coefficients that will be 
compared for each of these situations are P1Pl, P1P2, 
S1S1,SlP2. 

The situations that will be looked at are (a) 
isotropic medium over a transversely isotropic medium 
(b) transversely isotropic medium over an isotropic 
medium (c) transversely isotropic medium over a trans- 
versely isotropic medium (d) liquid medium over a 
transversely isotropic medium. The parameters of the 
media are tabulated below. (The velocities given 
should not be taken as absolute as it is only their 


ratios that are of significance.) 
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Description of the Media 
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PIP1 reflection coefficients 
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P1P2 transmission coefficients 
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S1S1L reflection coefficients 
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Figure A9 
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